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Abstract 

This paper treats the maximum queue length M , in terms of the number of 
customers present, in a busy cycle in the M/G/1 queue. The distribution 
of M depends both on the service time distribution and on the service 
discipline. Assume that the service times have a logconvex density and 
the discipline is Foreground Background (FB). The FB service discipline 
gives service to the customer(s) that have received the least amount of 
service so far. It is shown that under these assumptions the tail of M is 
bounded by an exponential tail. This bound is used to calculate the time 
to overflow of a buffer, both in stable and unstable queues. 

1 Introduction 

In a stochastic process often the extreme values rather than the usual values are 
of great interest. Large values of the queue length may ask for extraordinary 
measures such as the allocation of auxiliary storage space. In finite waiting room 
systems a natural question is: what is the probability that the queue length will 
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exceed a specific buffer size in a certain time period? Unlike the workload 
process, the queue length process is determined by the service discipline. 

In this paper we study the maximum queue length A/, measured in number 
of customers present in the system, in a busy period. The busy period maximum 
was studied by Cohen 5 , who gave an integral representation of the distribution 
of M in the M/G/1 FIFO queue. For the M/M/1 queue with load p < 1, relation 
(2.50) in Cohen 5 yields a simple expression for the exceedance probabilities: 



There is a growing interest in models with heavy-tailed service times, since 
statistical data analysis has provided convincing evidence of heavy-tailed traffic 
characteristics in high-speed communication networks. For heavy-tailed service 
times, the performance of the FIFO discipline decreases and it is natural to 
consider non-FIFO service disciplines. 

In this paper the service times have a log-convex density /, i.e. log / is con- 
vex. Both heavy-tailed distributions and light-tailed distributions may occur. 
Righter [Sj has shown that for queues with log-convex densities the Foreground 
Background (FB) service discipline minimises the queue length, see Theorem|21 
below. The FB service discipline gives service to the customers that have re- 
ceived the least amount of service so far. If there are n such customers, then 
they are served simultaneously at rate 1/n. When the age of a customer is 
the amount of service he has received, the FB discipline gives service to the 
'youngest' customer (s) present in the system. 

Our main result is that in the M/G/1 queue with a log-convex density oper- 
ating under the FB service discipline the tail of the maximum queue length M 
in a busy period is bounded by an exponential tail: 



where p is the load of the system. 

Interestingly, the upper bound does not depend on the precise form of the 
distribution. Note furthermore that (jJl is comparable to the precise value of 
the exceedance probabilities for exponential service times in J^l- 

By the regenerative structure of the queue length process, the maximum 
queue length over a busy period is related to the maximum over the time interval 
[0,t] for t ^ oo, see the survey article Asmussen Using the upper bound 

we show that in case of service times with heavy-tailed log-convex densities. 



P(M>n)=p"(l-p)/(l-p"+')- 
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P{M >n)<p", 
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the time to overflow of a buffer is of another order in the FB queue than in the 
FIFO queue. This illustrates the idea that using the FB discipline instead of 
FIFO may increase the performance of the queue considerably in case of heavy 
tails. 

The class of log-convex densities contains both heavy-tailed distributions 
like Pareto distributions (with density f{x) — (1 + x)^°' for a > 1) and WeibuU 
distributions with tail 1 — F{x) — exp(— a;'^) for f3 G (0, 1) and light-tailed dis- 
tributions like the gamma distribution with density f{x) = c{a, (3)x^ exp(— ax) 
for /3 > 0, a^c[a^l3) > 0. Equivalently, distributions with a log-convex density 
may be characterised by a decreasing likelihood ratio and have been called DLR 
distributions in the literature. 

Borst et al. [3] and Zwart and Boxma XI. discuss service disciplines in the 
case of heavy-tailed service times. For more results on the FB discipline, also 
known as FBPS or LAST, see Kleinrock [7j and the survey paper Yashkov jlOj . 

The paper is organised as follows. In Section [3 we establish the notation 
and prove the basic proposition. The inequality ^ is proved in Section |31 
and a sharper bound depending on the Laplace transform of the service-time 
distribution is established. In Section 01 we use (O to obtain asymptotics for 
the maximum queue length M{t) over the interval [0,<] for t — > oo and apply 
these to study the overflow time of the buffer in a stable queue. In SectionElwe 
describe a coupling that allows us to apply the results of Section 0] to the time 
to overflow in an unstable queue. 

2 Preliminaries 

In this section we establish a relation between the maximum queue length in 
the M/G/1 FB queue and in a related queue. This relation is the basis of the 
theorems in the next section. 

Consider an M/G/1 FB queue with Poisson arrival rate A and i.i.d. service 
times Bi, B2, ■ ■ ■ with distribution function F. Let B denote a generic service 
time and let the load p = XEB satisfy p < 1. Assume the queue is empty at 
time and let M be the maximum queue length in the first busy period. In 
order to determine an upper bound for the tail probabilities P{M > n), we 
compare the M/G/1 FB queue with a queue with the same arrivals and service 
times, but with an alternative service discipline, FB* say, defined as follows: 
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• the customer that starts the busy period has the lowest priority, and is 
only served when there are no other customers present in the queue 

• the other customers are served according to the FB discipline. 

Let M* denote the maximum queue length in a busy period in the M/G/1 FB* 
queue. The basic idea of this paper is contained in the following proposition 
that relates the systems FB and FB*. 

Proposition 1 For the discipline FB* above 

P{M* > n) = 1 - £;e-^^(^>"-i)-B. 

Proof By definition of the discipline FB* , the busy period consists of two (dis- 
joint) types of periods: periods in which only the customer that started the busy 
period is served and periods during which only other customers are served. The 
latter we call sub-busy periods. Let the random variable K denote the number 
of sub-busy periods in the busy period. The sub-busy periods are distributed as 
M/G/1 FB busy periods with one additional customer that is not served. Let the 
maximum queue lengths in the sub-busy periods be Mi + 1, M2 -|- 1, . . . , + 1. 
These are i.i.d. and have the same distribution as M -|- 1. 



1 Mi + 1 1 M2 + I 1 
I 1 1 1 1 1 

exp(A) exp(A) 

Figure 1: A realisation of the busy period in the M/G/1 FB* queue with K = 2. 

Conditional on the first service time Bi = x, the number of sub-busy periods 
_ftr is a Poisson distributed random variable with parameter Aa;. Indeed, K is 
the number of times the service of the first customer is interrupted. By the 
memoryless property of the arrival process, the times between two interruptions 
are independent and exp(A) distributed. 

Given that K = k and Bi = x the maximum queue length M* is the maximum 
of the k independent variables Mi + 1, . . . , M^ + 1 and hence 

P(M* <n\K = k,Bi=x) = {P{M < n))^ 

Since K is Poisson distributed for given x we find 



P(M* <n,K = k\Bi=x)= {P{M < n))'"^^^e 



-Ax 
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Summing over k yields a simple expression 



oo 



A:! 



P{M* <n \ Bi=x) {P{M < n)) 



e 



\xP{M>n-l) 



and 



P{M* <n)= e 
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The result now follows. 



□ 



3 Theorems 

In this section we derive exponential upper bounds for the tails of the maximum 
queue length M by combining Proposition ^ with Theorem |3 below. 

The following theorem is Theorem 13. D. 8 of Righter |H] in a form adapted to the 
present setting. It states that for service time distributions with a log-convex 
density, the FB discipline is optimal for minimising the queue length in a strong 
sense. 

As is usual, we restrict attention to service disciplines that cannot look into 
the future. This means that each instant the choice of the customer to be served 
is measurable with respect to the filtration {J-t,t > 0} where JFj is the tr-algebra 
generated by the arrival and departure times and ages of the customers up to 
time t. 

Theorem 2 (Righter) Let tt be a service discipline and let Nj:B{t) and iV^(i) 
denote the queue lengths at time t in G/GI/1 queues with disciplines FB and 
TT, respectively, under the same interarrival and service time distribution. If 
the service time distribution has a log-convex density, then there exist processes 
{NfB{t),t > 0} and {N.^{t),t > 0} such that {NFB{t),t > 0} = {A^fbW,* > 0}, 
{N^(t),t > 0} = {N^{t),t > 0} and 

PiNfBit) <N^{t),t> 0)^1. 

On applying Theorem|5]to Proposition^it is seen that P{AI > n) < P{M* > n). 
So we find 

Theorem 3 Let M be the maximum queue length in the busy period in the 
M/G/1 FB queue. IJ the service time distribution has a log-convex density, then 
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the exceedance probabilities r„ := P{M > n) satisfy 

r„+i < n = 0,l,2,.... (3) 

Theorem 4 Let M denote the maximum queue length in the busy period in 
the M/G/1 FB queue with workload p — XEB < I. If the service times have a 
log-convex density, then r„ = P{M > n) satisfies 

rn<p", n = 0,l,.... (4) 

In fact, ro = 1, ri = 1 — Ee^^^ and 

Tn < Prn~l, 71 > 1. (5) 

Proof By Theoreni|31we have r„+i < 4'{rn), where 0(a;) = 1 — i?exp(— ABx). 
Furthermore (/)(0) = and 

<l)'{x) = E\Beyi^{-\xB) <\EB = p <l. (6) 

Hence < 0(r„) = J^" (l)'{u)du < prn- This is lO, which imphes Q by 

induction since ro = 1. The probabihty ri may be computed exactly. The first 
interarrival time U is exp(A)-distributed and 

ri = P{M > 2) = P{Bi >U) = l-P{U > Bi) = 1-1 e-^''dF{x) = l-Ee^^^. 

Jo 

This proves the theorem. □ 

The approximation 1 — E ex.p{—XBx) w px in the proof of Theorem^ is quite 
good for X close to zero, since p is the slope of c/){x) = 1 — E exp{— XBx) in a; = 0. 

The following corollary to Theorem |2| gives even sharper upper bounds for 
P{M > n). 

Corollary 5 Let qo = I and set = 1 — Ee^^^^"^ for n > 1. Then 

P{M > n) <qn < {I- Ee-^^)p'^~\ n>l. 

Proof By Theorem^ induction and the fact that 1—E exp{—XBx) is increasing 
in x, we have for n > 

P{M >n)<l- £;g-AP(M>n-i)B < ^ _ Ecxp{-Xqn^iB) = g„. 
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Inequality © and induction yield qn+i ^ PQn 

Example Numerical calculations show that the bounds qn in Corollary may 
be significantly better than the bounds in Theorem 01 For service times with 
the Pareto density 3(1 + x)~^ and arrival rate A — 1.8, the load is p = 0.9 and 
(l-£;exp(-AB))p97qioo«7.5. 

If the service times in one queue are stochastically smaller than those in an- 
other queue with the same arrival rate and service discipline, the maximum 
queue length in the first queue is stochastically smaller than that in the second 
queue. Hence the upper bounds in this section also hold for service times that 
are stochastically smaller than a service time with a log-convex density. The 
next corollary states this idea for the Pareto distribution. 

Corollary 6 Consider an M/G/1 FB queue with arrival rate A and generic 
service time B. Let M denote the maximum queue length in a busy period. Let 
a > 2 and c > 0. // 

P{B > x) < x>0, 
^ ^ ~ (1 -f ex)" 

then P{M > n) < 6*" for all n > where 9 = A/((a - l)(a - 2)c'^). 



4 Asymptotics for the maximum queue length 
over an interval 

In this section we present asymptotics for M{t), the maximum queue length 
in the interval (0, i), for t — > oo. These are applied to calculate the time to 
overflow in a finite-buffer system. 



A stochastic process X{t) is called a regenerative process if there exists a (pos- 
sibly delayed) renewal process with epochs < Tq < Ti < ■ ■ ■ such that the 
cycles 

{X{t + Ti_i)}o<t<Ti-Ti_i 

are i.i.d. for i = 1, 2, . . .. 
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Consider an M/G/1 FB queue with arrival rate A, i.i.d. service times with log- 
convex density /, generic service time B and workload p < 1. The initial state is 
arbitrary. Since the interarrival times are memoryless, the queue length process 
is a regenerative process, where a cycle consists of an idle and a busy period. 
The expected cycle length fi is given by ^ = A^^ + EB/{1 — p) = A^^(l — p)~^ ■ 
Let M{t) denote the queue length over the time interval [0, t\. Proposition VI. 4. 7 
in Asmussen 2 then states that for t — > oo 

sup \P{M{t) <x)~ P{M < x)*/^! -> 0. (7) 

X 

Example (Buffer overflow) We are interested in the time to overflow of a 
buffer of size d > 1. Denote by td.p = mi{t : P{M{t) > d) > p} the lower pth 
quantile for the (random) time to overflow. Since td,p ^ cxd as d ^ oo, we have 
by 10 that 

P{M{td^p) <d)= P{M < dY^'■-/^'{l + o(l)) 

as d ^ oo. Theorem 0] then yields 

/^log(l -p) ^log(l -p) 
^ \ogP{M <d)- log(l - pd) ^''> 

where the asymptotic equality holds as o? ^ cxd. The RHS of (jH)) may be sharp- 
ened by using q„ instead of p", where q„ is defined in Corollary El 

We compare the time to overflow in the FB and the FIFO queue. Set p = ^. 
Then td^p is the median of the time to overflow. Suppose the service times have 
a Pareto distribution with distribution function 1 — {x + 1)~^ and the arrival 
rate A is 1.8, so that the load of this queue is p = 0.9. Time is measured in mil- 
liseconds. Let the buffer size d be 1000. Using the asymptotic inequality 0, we 
flnd for the FB discipline a value of td^p larger than 10"*^. This is approximately 
10^^ years. 

Now consider the same queue with the FIFO discipline. Occasionally cus- 
tomers with very large service times arrive. Such customers may cause a buffer 
overflow. For large values of d, the probability that a customer with service time 
at least d/A + Su, with a ~ ^ dj A, will cause a buffer overflow is larger than 
0.99. The probability that such a customer arrives during the interval [0,t] is 

p(t) = 1 - exp(-At(l - F(d/A + 3VrfA)))- 

Solving p(ii) = 1/2 we find t\ < 10^, approximately one day. 
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For the M/G/1 FB queue with load p = 0.9 and WeibuU service time distri- 
bution F{x) = 1 — exp(— a;''), the value of to is larger than 10'*^ for (3—1/4: and 
larger than 10'*^ for /3 = 1/2. For the FIFO queue with the same service time 
distribution and load, ti is smaller than lO'^ for (3 = 1/4 and approximately 10^° 
for P = 1/2. 

5 Buffer overflow in unstable queues 

In this section we apply the results of section 4 to study the time to overflow of 
the buffer in unstable queues. 

When the condition p < 1 is violated, the workload asymptotically grows 
with rate p — 1. Balkema and Verwijmeren showed that in the overloaded 
M/G/1 queue under the FIFO discipline the queue length asymptotically grows 
linearly with rate A, where A is the arrival intensity. The queue length in the FB 
queue grows linearly as well, but with a smaller rate. By the priority rule of the 
FB discipline, customers with service time less than the critical value c*, where 
c* = inf{c : XE{BAc) > 1}, will a.s. leave the queue, since they are not hindered 
by customers with long service times. Customers with service time larger than 
c* have a positive probability of being stuck in the queue forever. Here we are 
interested in the time to overflow of a buffer of size d in an overloaded queue. 

Consider a queue with heavy-tailed service times and a small arrival rate. 
Assume the stability condition p < 1 is violated. This happens because of 
the mass in the far right tail of the service time distribution. By modifying 
the right tail, we shall obtain an alternative service-time distribution that does 
satisfy the condition p < 1. The maximum queue length in the alternative 
queue is comparable to that in the original queue for very long periods of time 
if the arrival process is light. We shall construct the alternative service times in 
a way such that they have a log-convex density. Then Theorem 0] and O are 
be applied to bound the tail of the maximum queue size. 

For the comparison between the two queues we use a coupling argument 
described in Subsection 15.11 below. Subsection 15.21 illustrates these ideas with 
an example. 

Extreme value behaviour of M in the M/G/1 FIFO queue was found by 
Cohen under the assumption of the existence of an exponential moment. 
Heavy-traffic limits in this case were considered by Iglehart [Hj. Serfozo jH] 
studied the behaviour of M{t) both in stable and unstable GI/G/1 FIFO queues. 
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5.1 The coupling argument 

Theorem 7 Let Mp{t) and MQ{t) he the maximum queue lengths over the 
interval [0, i] in two M/G/1 FB queues with service-time distribution functions 
F and G, and the same arrival rate A. Assume FAp — G/\p for somep G (0, 1). 
Then 

Mpit) <st Mait) + Kp(t) 
where Kp{t) is a Poisson process with rate A(l —p). 

Proof Let Np{t) and Nc{t) be the queue lengths at time t. We prove that 

Npit) <st Nait) + Kpit) (9) 

where Kp{t) is a Poisson process with rate A(l — p). Using the fact that Kp{t) 
is non-decreasing, the theorem then follows from 

Mpit) = nvAyi{NF{s)} <st max{7VG(s) + Kp{s)} < Mait) + Kp{t). 

s<t s<t 

The coupling is standard: let J/i, i72, ■ • • be a sequence of i.i.d. random variables 
distributed uniformly on [0, 1]. Define the sequences Bi, B2, ■ ■ ■ and Bj}" , i?J , . . . 
by setting Bk = F-\Uk) and B* = G-^{Uk) for fc = 1, 2, . . .. Then Bk and B* 
have distribution functions F and G respectively, for all k, and 

BkAp = BtAp, fc = l,2,.... 

Now let Q and Q* be the M/G/1 FB queues with the same arrival process 
with rate A and service times introduced above. Let N{t) and N*{t) be the 
queue lengths at time t and let Np{t) and Np{t) be the part of the queue 
lengths formed by customers younger than F~^{p). By the coupling of the 
service times and the fact that the FB discipline favours customers younger 
than F^^{p) over customers older than F~^{p), we have Np{t) — N*{t). Let 
Kp{t) be the Poisson arrival process of customers with service time larger than 
F-^{p). Then N{t) < Np{t) + Kp{t) and N*{t) < N*{t). Hence 

Npit) = N(t) < N*it) + Kp{t) = Nait) + Kp{t). 

We conclude the proof by observing that the rate of Kp{t) is A(l — p). □ 

Note that the random variables Mc{t) and Kp{t) in Theorem 13 are dependent. 
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5.2 Example 

In this subsection we compute an upper bound for the time to overflow of an 
M/G/1 FB queue with Pareto service density f{x) = l/(a; + 1)^, so that p = oo. 
For a > 1 define 



Then f{x) — ga{x) for x < a and it may be shown that ga is a probabihty 
density which is continuous and log-convex. Now let M{t) and Ma{t) denote 
the maximum queue lengths over the interval [0, i] of two M/G/1 FB queues 
with service-time densities / and ga respectively, and with the same arrival rate 
A. Setting p = 1 — (a + 1)^^, we have by Thcorcm[3 



where Ka{t) is a Poisson process with rate A/(l + a). Hence for any xi, a;2 > 
such that xi -\- X2 ^ X we have 



Setting a = 10 and A = 0.01, we find for the queue with density ga that 
p = XEB = A(log(a + 1) + 1) = 0.9 . . .. Since for this queue p < I and ga is 
log-convex, we may apply lO and lO to Ma{t). For a buffer of size d — 1000 
and p — 0.01 we find that td~i,p — 3 • 10^^(1 -I- e), where the error e = s{d,p) is 
determined by the asymptotic inequality in 0. Setting xi = d — 1 and X2 ~ 1 
in (|11() then yields that the probability of a buffer overflow in [0, t] with t = 10'*° 
is smaller than 0.02. 

We conclude that for unstable M/G/1 FB queues, the time to overflow may 
be very large, provided that the arrival rate is low. 
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M{t) <st Ma{t) + Ka{t) 



(10) 



P{M{t) >X)< P{Ma{t) > Xi) + P{Ka{t) > X2). 



(11) 
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